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Abstract 

We study effects of charged Higgs boson exchange in B —* DtT> t . The Yukawa 
couplings of Model II of two-Higgs-doublet model, which has the same Yukawa 
couplings as MSSM, is considered. We evaluate the decay rate including next-to- 
leading QCD corrections and estimate uncertainties in the theoretical calculation. 
Our analysis will contribute to probe an extended Higgs sector at B factory exper- 
iments. 
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1 Introduction 



The minimal supersymmetric standard model (MSSM) |T| is one of the most attractive 
models beyond the standard model (SM). In the MSSM, two Higgs doublets are introduced 
in order to cancel the anomaly and to give the fermion masses. The introduction of the 
second Higgs doublet inevitably means that a charged Higgs boson is in the physical 
spectra. So, it is quite important to study effects of the charged Higgs boson. 

Here, we study effects of the charged Higgs boson on the exclusive semi-tauonic B 
decay, B — ► Dtu t , in the MSSM. In our previous works @, [3|, we calculated the decay 
rate of B — > Dtv t including the effect of charged Higgs boson exchange in the leading 
logarithmic approximation and the heavy quark limit. In this work, we show the decay rate 
with QCD corrections up to the next-to-leading order (NLO). The NLO corrections are 
necessary to estimate theoretical uncertainties coming from short- distance calculations in 
the ratio of the decay rates (see below). In addition, these corrections may cause dominant 
uncertainties for the q 2 distribution |4[] and the r polarization [[|. 

In a two-Higgs-doublet model, the couplings of charged Higgs bosons to quarks and 
leptons are given by 

C H = (2v / 2G F ) 1/2 [Xu L V KM M d d R + Yu R M u V KM d L + ZV L M[l R ] H + 

+h.c. , (1) 

where M u , Md and Mi are diagonal quark and lepton mass matrices, and Vkm is Kobayashi- 
Maskawa matrix ||. In the MSSM, we obtain 

X = Z = tan p , Y = cot (3 , (2) 

where tan/3 = V2/V1 is the ratio of the vacuum expectation values of the Higgs bosons. 
Since the Yukawa couplings of the MSSM are the same as those of the so-called Model II 
of two-Higgs-doublet models ||, the above equations and the following results apply to 
the latter as well.Q 

With these couplings, it turns out that the amplitude of charged Higgs exchange in 

B — > Dtv t has a term proportional to m^tan 2 (3. Therefore, the effect of the charged 

Higgs boson is more significant for larger tan (3. 

1 It is known that SUSY loop effects in Eq.(fJ) arc significant for large tan/3 0. The dominant 
effect in the b — > crv T decay || comes from the SUSY-QCD correction to the bottom quark mass ||. 
Once this correction is taken into account, the Yukawa couplings of the MSSM is no longer the same as 
Model II of two-Higgs-doublet models. This effect cannot be ignored in order to study the MSSM Higgs 
sector. However, we omit it in this work because our aim is to clarify low-energy QCD corrections and 
uncertainties, which are universal and model-independent. The SUSY loop effects on B — > DtD t will be 
discussed elsewhere 



1 



Formula of the decay rate is described in Sec. 2. In Sec. 3, we give hadronic form 
factors including next-to-leading QCD corrections. In Sec. 4, we show our numerical 
results. Our conclusion is given in Sec. 5. 



2 Formula of the decay rate 

Using the above Lagrangian in Eq. (|l|) and the standard charged current Lagrangian, we 
calculate the amplitudes of charged Higgs exchange and W boson exchange in B — ► Drv r . 



The W boson exchange amplitude is given by JTT 

G p 

71 



M X /(q 2 ,x) w = ^V cb J2vx w Ltm w , (3) 



where q 2 is the invariant mass squared of the leptonic system, and x = p B • p T /m 2 B . The 
r helicity and the virtual W helicity are denoted by A r = ± and \w = ±, 0, s, and 
the metric factor r]\ w is given by rj± = rjo = —rj s = 1. The hadronic amplitude which 
describes B — > D W* and the leptonic amplitude which describes W* — > rv are given by 

H{ w {q 2 ) = e;(\ w ){D(p D )\ry%\B(p B )) , (4) 
L x x l(q 2 ,x) = e,(\ w )(T(p T ,\ T )v T (p u )\fY(l-l 5 )v T \0) , (5) 

where e^(Xw) is the polarization vector of the virtual W boson. 
The charged Higgs exchange amplitude is given by || 



M^(q 2 ,x) H = ^V cb L^ 



M 2 H R M 2 H L 



(6) 



Here, the hadronic and leptonic amplitudes are defined by 



H s R , L {q 2 ) = (D(p D )\c(l± l5 )b\B(p B )), (7) 
L^(q 2 ,x) = (t(p t ,\ t )9 t (p„)\T(1-j 5 )v t \0) . (8) 

This leptonic amplitude is related to the W exchange amplitude as 

\fcP- 
m T 

Details of the hadronic amplitudes for the W exchange and the charged Higgs exchange 
are discussed in the next section. 
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Using the amplitudes of Eqs. (|3]) and (^|), the differential decay rate is given by 

dT ^2| T /J2^4 



dq 2 



128vr 3 m| L\3 y 3 V V o; 
+m^(m 6 tan 2 /? + m c )^^-^) 2 ] , (10) 



where Q± = ± m^) 2 — g 2 and = x/T— m 2 T jq 2 . Note that if tan /3 > 1, in which we 
are interested, this decay rate is practically a function of tan/3/M# because the second 
term in the coefficient of H S R is negligible for m^tan 2 f3 ^> m c . 

3 Hadronic form factors including QCD corrections 

Here, we evaluate the hadronic amplitudes in Eq. (|j) and Eq. (0) in order to obtain the 
decay rate numerically. These amplitudes are given in terms of hadronic form factors: 

(D(p D )\cYb\B(p B )) = ^^B[h + {w){v + v'Y + h_{w){v - v'Y] , (11) 
(D(p D )\cb\B(p B )) = sM^(l + w)h s (w) , (12) 

where v = Pb/i^b, v' = Pd I m D an d w = v ■ v' = (m 2 B + m 2 D — q 2 )/2rriBm D |12| . 

In the heavy quark limit and in the leading logarithmic approximation, these form 
factors h±{w) and h s (w) are given as 

h+(w) = h s (w) = £(w) , /i_H=0, (13) 

where £ (w) is the universal form factor ||TB| . 

Now, we consider QCD corrections beyond LLA and calculate these form factors up 
to the next-to-leading order. Then, these form factors are given as 

h+(w) = {c 1 (w)-(^){C 2 (w) + C 3 (w))}t(w), (14) 

h-H = -{^{C 2 (w)-C 3 (w))}aw), (15) 
h s (w) = C s (w)£(w). (16) 



Explicit formula of coefficients for W exchange, Cj (i = 1 <~ 3), are given by Neubert [14 
The coefficient for charged Higgs exchange, C s , is given as 

C s (w) = A(w)C s {w) , (17) 
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where 

A(w) 
C s {w) 



1 + 



a s (m b ) 

a s (m b ) - a s (m c 



[Z + 2) + 



7T 



a s [m c \ 



7T 



Z(«;) + 2 + -(/H-rH) 



2a s (m) 
H ~ 9s{z,w) . 



3tt 



(18) 



(19) 



m is some average mass of m& and m c and other functions and constants in Eq. 
and Eq. (|1^) are given in Ref.fll4l and Appendix. We have used the MS scheme in our 
calculations. 



The form of £ (w) is strongly constrained by the dispersion relations as [|T5 
£0) ~ 1 - 8pjz + (51.pl - 10.)z 2 - (252.pl - 84.)z 3 , 



(20) 



where z = (\Jw + 1 — V2)/(y/w + 1 + y/2). We obtain the slope parameter pi as 

pi = 1.33 ± 0.22 , (21) 



from the experimental data of B — > D*ev e [16 



4 Numerical results 

We consider the ratio of decay rates, 

B = T{B - D ™ t) (22) 
T(B -> D/i^)5M ' 1 J 

where the denominator is the decay rate of B — > Dpv^ in the SM. Uncertainties due to 
the form factors and other parameters tend to reduce or vanish by taking the ratio. 

Fig. 1(a) is the plot of our prediction of the ratio in Eq. ( p2"D as a function of R, which 
is defined by R = m^y tan/3/m#. Here, we do not show the error in the slope parameter 
pi and we take pi = 1.33. The dashed lines show the MSSM and SM predictions without 
QCD corrections. The lines with narrow shaded regions show the predictions including 
QCD corrections with Aj^g being varied between 0.15 GeV and 0.25 GeV. 

In Fig. 1(b) we show the ratio of B in the MSSM with and without QCD corrections, 
5 (with QCD corrections)/!? (without QCD corrections), as a function of R. The solid 
line is the ratio with = 0.25 GeV and the dashed line is that with A^jg = 0.15 GeV. 

From Fig. 1, we expect that theoretical uncertainties from higher order QCD cor- 
rections are at most a few percents in the ratio of decay rates and, as seen later, the 
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ratio 




Figure 1: (a) The ratio B as a function of R at p\ = 1.33 in the MSSM and the SM. The 
lines with shaded region are obtained by using Aj^g = 0.15 ~ 0.25 GeV and the dashed 
lines show the predictions without QCD corrections, (b) The ratio of B with and without 
QCD corrections, £?(with QCD corrections)/!? (without QCD corrections), as a function 
of R. The solid and the dashed lines are the ratios with Atts- = 0.25 and 0.15 GeV. 



theoretical uncertainties from QCD corrections are much smaller than those from the 
error of p\. 

Fig. 2(a) shows our prediction of the ratio B with QCD corrections as a function of 
R. Here, we take the error in the slope parameter into account and use Aj^g- = 0.25 
GeV. The shaded regions show the MSSM and SM predictions with the error in the slope 
parameter p\ in Eq. (j2"T|) . As mentioned before, from Fig. 2(a) and Fig. 1, we see that the 
theoretical uncertainty from the error in the slope parameter is dominant over that from 
QCD corrections. As seen in Fig. 1(a), when R is about 35, the ratio B in the MSSM is 
the same as the one in the SM. 

In Fig. 2(b), we also show the ratio, 

B= m^D T! > T ) 

Y{B -> Dfiu^sM 

with QCD corrections, similar as Fig. 2(a), but its denominator is T(B — > Dpu^sM, 
which is integrated over the same q 2 region as the r mode, i.e., m 2 < q 2 < (ms — wid) 2 - 
From Fig. 2(b), we observe that the ratio B has less theoretical uncertainty and we expect 
a better sensitivity than the ratio B in Fig. 2(a). 

Now, we consider a decay distribution defined by 

d(w) = (w 2 - 1 — JL _ T)/ , (24) 

T(B -> Dpv^sM/dw 
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Figure 2: The ratios B and B with QCD corrections as functions of R in the MSSM and 
the SM. The shaded regions show the predictions with the error in the slope parameter p\ 
in Eq. (pl|) . The flat bands show the SM predictions, (a) B: the decay rate normalized 
to T(B — > Dfiu^sM- (b) B: the same as (a) except that the denominator is integrated 
over m 2 < q 2 < (mg — mo) 2 - 



where w = (m 2 B + m 2 D — q 2 ) /2mBrri£>. 

In Fig. 3, we show d(w) in the SM and the MSSM with different values of R. The 
dashed lines show the MSSM and SM predictions without QCD corrections. The lines 
with shaded regions show the predictions with QCD corrections and the shaded regions 
are given by using Ajyjg = 0.15 ~ 0.25 GeV. The theoretical uncertainty from the error 
in the slope parameter is canceled out in this quantity. Thus, QCD corrections become 
dominant uncertainties in the theoretical calculation. From Fig. 2, the ratio B in the 
MSSM becomes the same as the one in the SM when R ~ 35. But, in Fig. 3, we find that 
the behavior of d(w) in the MSSM with R = 35 is considerably different from that in the 
SM. Therefore, we can distinguish the SM from the MSSM by investigating w distribution 
even if R ~ 35. 



5 Conclusion 

As seen in our numerical results, the branching ratio of B — > Drv T is a sensitive probe of 
the MSSM-like Higgs sector. So, if B — > Dtv t is observed at a B factory experiment, a 
significant region of the parameter space of the MSSM Higgs sector will be proved. This 



is complementary at the Higgs search at LHC [[17] . 

In the branching ratio, the theoretical uncertainty from QCD correction is much 
smaller than that from the error in the slope parameter p\. However, in w distribu- 
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Figure 3: The w distribution d(w) in the SM and the MSSM with different values of R. 
The lines with shaded regions are given by using Aj^g = 0.15 ~ 0.25 GeV and the dashed 
lines show the prediction without QCD corrections. 



tion, the theoretical uncertainty from the error in the slope parameter p\ is canceled out 
and, therefore, it is important to consider QCD corrections. 

In this work, we have not taken 1/m corrections into account. They may be as 
significant as QCD corrections in the w distribution. However we expect that their effects 



are smaller than the uncertainty from the error of p\ in the branching ratio |10 



Finally, the r polarization in B —> Drv r is also expected to be a good probe of charged 
Higgs boson. The theoretical uncertainty from the error in the slope parameter becomes 
very small in the r polarization ||. QCD and 1/m corrections to this quantity will be 
addressed elsewhere. 
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Appendix 



Functions in Eq. fll8|) and Eq. ( fL9|) are given by 

1 



r(w) 

/H 

a L (w) 
Z{w) 
Z 



w- 



ln(u> + Vw 2 — 1) 



wr(w) — 2 



L 2 (l-wi)-L 2 (l-w1 



1 — + z 2 
wr(w) — 1 



2Vw 2 - 1 L 
L 2 (l — — L 2 (l — zw 

[w 2 — l)r(w) + (w — z) In z 



27 
4 

~9 
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L54 + 12 + 9 + ^ W '\ l wr ( w ) Q^ w > ' 

9403 



-7T 2 - 



225 7500 ' 
where z = m c jm^ w± = w ± Vw 2 — 1, 



ln(l -t) 



dt 

o * 



is the dilogarithm and, 
I(w) 



dip 



ip coth xp — 1 



x j^coth 2 <p + 



sinh </? cosh ip sinh 



In 



sinh 2 ip - sinh 2 r/> ~ sinh ip . 
in terms of the hyperbolic angle <p defined by w = coshy?. 



(A.l) 
(A.2) 

(A.3) 
(A.4) 
(A.5) 
(A.6) 



(A.7) 



(A.8) 
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